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Abstract 

Parameter-dependent statistical properties of spectra of totally connected irregular quantum 
graphs with Neumann boundary conditions are studied. The autocorrelation functions of level 
velocities c(x) and c(cj,x) as well as the distributions of level curvatures and avoided crossing 
gaps are calculated. The numerical results are compared with the predictions of Random Matrix 
Theory (RMT) for Gaussian Orthogonal Ensemble (GOE) and for coupled GOE matrices. The 
application of coupled GOE matrices was justified by studying localization phenomena in graphs' 
wave functions *$>(x) using the Inverse Participation Ratio (IPR) and the amplitude distribution 
P(*(a?)). 
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Hamiltonians H(X) of many chaotic systems depend on external parameters X. Depenc - 
ing on the system there is a big variety of parameters X. It may be the potential energy [l], 



the magnetic field Q, 3, 0, J Q, H, |8[, 
lengths of the bonds of quantum graphs 



the shape of billiard's boundary 



10| or the 



11]. For parameter dependent systems the energy 
levels are determined through the eigenvalue problem H(X)4>i(X) = E i (X)4> i (X). It is 
conjectured that for quantum chaotic systems, when the variation of the parameter X does 
not change the symmetry of the system, statistical properties of the level dynamics should 
be universal {2], I3], O, O, Q, Q] • In order to investigate universalities in the level dynam- 
ics several parametric statistics can be used, for example, autocorrelation functions of the 

n * n 

level velocities [5j , the level curvature distribution and the avoided crossing distribution |14J | . 
Parametric dynamics of random matrices was also investigated la. 19. 171 18, 19, 

Quantum graphs have attracted much attention in recent years |11| . 1211 . 1221 . |2j 
A lot of properties of quantum graphs have already been studied. For example, spectral prop- 
erties of graphs were studied in the series of papers by Kottos and Smilansky [271], where the 
authors showed that quantum graphs are excellent paradigms of quantum chaos. The auto- 
correlation functions of level velocities for graphs with five vertices with and without time 




reversal symmetry were studied in the paper 



28] . However, many other properties of graphs 



require much thorough consideration. In this paper we study the autocorrelation functions 
of level velocities as well as the distributions of level curvatures and avoided crossing gaps 
for quantum graphs. 

Graphs can be considered as idealizations of physical networks in the limit where the 
widths of the wires are much smaller than their lengths. Among the systems modelle d by 



graphs one can find, ag., electromagnetic optical waveguides |2j3|, |3C| , quantum wires 
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mesoscopic systems 



33 



341 ]. or microwave networks 3a. 136] 



32|, 



A quantum graph (network) consists of n vertices connected by bonds. On each bond of 
the graph one-dimensional Schrodinger equation is defined (we assume that H = 2m = 1): 

d? 

(1) 



where k is the wave vector and the subscripts i, j denote the bond which connects two 
vertices with the numbers % and j. More detailed description of quantum graphs can be 
found in the paper jll| . 

In order to investigate the autocorrelation functions of level velocities one should unfold 
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both the energies E{ = kf and the parameter X. Energy levels of quantum graphs are 
unfolded by using the mapping 

£ = N^Et), (2) 
where N av (E) is the average number of states, which assures that the mean leve 



equal to unity. The average number of states N av (E) is given by the formula 27] 



spacing is 



AU^) = — + ^ (3) 

7T I 

where L is the total length of the graph. 

The parameter X was unfolded using the generalized conductance Cq(X) 

<W0-((W). 
where (. . .) means the average over the energy levels. The unfolded parameter is given by 



the following formula 371 ] : 



x 



C (X)dX, (5) 



x = 

where [X in ,X] is the interval of integration. After unfolding of the parameter X statistical 



properties of energy levels £»(x) of chaotic systems should be universal j^, 37j. Particulary, 
for quantum chaotic systems the parametric level velocities distribution Vj = d^/dx should 
be described by a Gaussian distribution jf], Q • 

The autocorrelation function of level velocities c(x) is defined as follows: 

c(x) = (lf(x)l?(x + x)Y (6) 



dx dx 

where the average is performed over the parameter x and over all the levels. This function 
is the measure of correlation of level velocities which belong to the same energy level. 
Another autocorrelation function investigated in this paper is c(u>, x) jlj] 

~, x) = S M ^te(x)-^(x + x)-^)f(x)f (x + x)) 
1 ' } E^(x)-£;(x + x)-u;)> ' [ ' } 

where the average is performed over the parameter x only. This function measures the 

correlation of level velocities for energy levels separated by a distance x in space and by a 

distance u in the unfolded energy. 

We study the autocorrelation functions of the level velocities c(x) and c(u, x) for quantum 

graphs. We consider only totally connected quantum graphs (every vertex is connected with 



the others) with no loops and no multiple bonds (every two vertices are connected by one 
bond only). The number of vertices n defines a graph size. For graphs considered in this work 
it was varied between n = 4 and n = 30. We impose the Neumann boundary conditions on 
each graph's vertex which imply continuity of the wave function and the probability current 
conservation on the vertices, n vertices of a totally connected quantum graph are connected 
by B = n[n — l)/2 bonds which gives B = 6 bonds for the graph with n = 4 vertices and 
B = 435 bonds for the graph with n = 30 vertices. 

The change of the bonds lengths of a graph was chosen to be an external parameter X. 
For the graphs with the even number of bonds the lengths of all the bonds were changed, 
while for the graphs with the odd number of bonds we changed the lengths of all the bonds 
except the arbitrary chosen one. The lengths of one half of the arbitrary chosen bonds of 
a graph were increased Lij(X) = Ly + X while the lengths of the other half of the bonds 
were decreased Li/ t f(X) = Li/ji — X. In this way the total length of the graph L = J2i<j Lij 
and the mean density of states were independent on the parameter X. The lengths of the 
graph's bonds were changed in 100 equally spaced steps dX: X = mdX (m = 1, 100). 

The lengths of the bonds L^j of each graph were chosen according to the formula 

L' 

Lij = ^ J-, > ( 8 ) 

where were uniformly distributed on (0, 1]. The total length of the graph was equal to 
one L — 1. For every size of a graph we realized 99 graph's configurations with different 
bonds lengths. To ensure the bonds lengths to be much larger than the wavelength A = 27r/ k, 
we imposed the condition L^(X) > 10A to be satisfied for all graphs configurations. 

For each configuration of the graph about 80 eigenvalues with the level numbers between 
822902 and 822985 were calculated. 

Figure [T] shows a typical behavior of unfolded energy levels 822930 < £j < 822960 as a 
function of the external parameter x for graph with n = 30 vertices. The parameter x was 
unfolded using the formula ([5]) . In Figure [T] one can see very complicated level dynamics 
with many avoided crossings. 

Parametric level velocities of quantum graphs were calculated using the finite difference 
method 

_ <9& &(x + rfx) -&(x) 
Vi ~ dx - dx • [J) 
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FIG. 1: Parametric energy level dynamics for totally connected quantum graph with n = 30 
vertices. 



The level velocities were rescaled using the variance cr v 



{{d^/dxf) 



(10) 



Figure [2] shows the velocity distribution P(v) for graphs with n — 6, 10 and 30 vertices. 
The results are averaged over 99 graphs configurations (approximately 8000 data points for 
each configuration). In all the cases presented in Figure [2] the velocity distribution P(v) 
calculated for quantum graphs is in good agreement with the Gaussian distribution (solid 
line), as it is expected for quantum chaotic systems. For the other graphs which are not 
presented in Figure [2] with the number of vertices between n = 6 and n = 30 good agreement 
with the Gaussian distribution is also observed. However, for the graphs with n = 4 and 
n = 5 vertices situation is different (results are not shown). While for the graphs with 
n = 5 the agreement with the Gaussian distribution is still quite satisfactory, for the graphs 
with n = 4 the velocity distribution is lower than the Gaussian curve for small values of 
level velocities. This result is not surprising, the graph with n = 4 is the simplest possible 
non-trivial fully connected graph one can construct and for some spectral statistics it could 
be too simple to exhibit behavior that is typical for quantum graphs. But when the number 
of vertices of a graph increases the agreement with the Gaussian curve improves. For this 
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FIG. 2: The parametric velocity distributions P(v) for quantum graphs with n = 6 (open circles), 
n = 10 (open triangles) and n = 30 (full circles) vertices, respectively. The numerical results are 
compared to the Gaussian distribution (solid line). 



reason in our further investigation we will concentrate mainly on graphs with n > 6 vertices. 

Figure [3] shows the autocorrelation function c(x) calculated for quantum graphs with 
n = 4, 5, 6, 7, 10, 20 and 30 vertices. In all cases the autocorrelation function c(x) was 
obtained by averaging over 99 graphs configurations. Numerical results are compared with 
the results of RMT for GOE (solid line). The numerical simulations of the parameter- 
dependent autocorrelation function c(x) within RMT for GOE were made using the following 



Hamiltonian model 10] 



H(X) = #isin(X) +H 2 cos{X), (11) 

where Hi, H 2 are two 500 x 500 matrices that are members of GOE. In our calculations 
the parameter X was chosen in 1001 equally spaced points in the interval [0, tt/8}. For each 
value of the parameter X we ran 99 realizations of Hi and H 2 . The eigenvalues obtained 
from the diagonalization of the Hamiltonian H(X) were unfolded by using the integrated 
average eigenvalue density for GOE matrices [3^]. The parameter X was unfolded using the 
formulas (4) and (5). 

As it is shown in Figure [3] the autocorrelation function of level velocities c(x) calculated 
for the graphs with n = 5 deviates in the downward direction from the RMT predictions for 



J , L 
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FIG. 3: The level velocity autocorrelation function c(x) for quantum graphs with n = 4 (crosses), 
n = 5 (full diamonds), n = 6 (open circles), n = 7 (full squares), n = 10 (open triangles), n = 20 
(open diamonds) and n = 30 (full circles) compared to the predictions of RMT for GOE (solid 
line). 



GOE (solid line) for most of the values of the parameter x, while for the graphs with n> 7 
the autocorrelation c(x) deviates in the upward direction from the RMT. The larger graph 
is, the bigger is deviation from RMT. Only for the graphs with n = 6 the agreement with 
RMT is quite good. For the graphs with n = 4 a strong deviation from RMT is observed for 
0.75 < x < 1.5, what is the evidence of a non-universal behavior of the system. It is worth 
mentioning that the difference between the results obtained for big graphs with n = 20 and 
n = 30 is very small what might suggest that for graphs with n > 30 the autocorrelation 
function c(x) converges to the limiting curve which is not given by the RMT predictions for 
GOE. 

The reason of the disa gree ment of c(x) from the RMT predictions is not exactly known. 
However, one can assume [39I- l4ol| that it is connected with non-ergodic structures of graphs' 
wave functions. In order to check this assumption we calculated the Inverse Participation 
Ratio (IPR) that is a measure of wave functions localization. 
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FIG. 4: The inverse participation ratio (IPR) calculated for quantum graphs (full circles) compared 
to the results of random wave hypothesis (solid line). The dotted line shows the minimum values 
{Imin = 1/2-B) °f IPR- In the inset we show the distributions P(^>(x)) of the wave functions ^f(x) 
for quantum graphs with n = 4 (crosses), n = 5 (full diamonds), n = 6 (open circles), n = 7 (full 
squares), n = 10 (open triangles), n = 20 (open diamonds) and n = 30 (full circles) compared to 
the Gaussian distribution (solid line). Distributions are presented in double logarithmic scale. 



In each bond of a graph a wave function can be written as [4 



(12) 



Then the inverse participation ratio for the N-th energy level can be defined in the following 

way: 

In = £2B ' i\ ' 2 - (13) 

l^2B \ a i,j I 

It takes the values between l m in = 1/25 for states which occupy each directed bond with 
the same probability and l ma x — 1 for a state which is restricted to a single bond only 
i.e., the greatest possible degree of localization [39J. Random Wave Hypothesis (RWH) [40 ] 
assumes Gaussian random fluctuations of the complex coefficients dij. Thus (\a.ij\ A ) = 2 
and {\aij\ 2 ) = 1. In this case Irwh ~ 1/B. In FigureHJwe compare the inverse participation 
ratio I n = (In) calculated for quantum graphs of different size (circles) with the RWH one 
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(solid line). The averaging in the last formula was performed over the energy levels and over 
graphs configurations. The dotted line in Figure H] shows the minimum value (I m i n — 1/2-E?) 
of IPR. Figure H] shows that for graphs with n > 5 the obtained results are bigger then RWH 
predictions Irwh- It means that wave functions of quantum graphs are less ergodic than 
the ones predicted by Random Wave Model. Furthermore, it is important to noting that 
the result obtained for the graphs with n = 5 and n = 6 vertices are the closest ones to the 
Irwh predictions. 

We also calculated the distributions P(\l/(x)) of the values of the graph wave function 
ty(x) 41j. The distributions P(\l/(x)) for graphs with n = 4,5, 6, 10, 20 and 30, in the double 
logarithmic scale, are presented in the inset in the Figure HI The distribution P(\l/(x)) is 
symmetric, therefore, we considered only positive values of the wave function ty(x) > 0. The 
inset shows that the distribution P(^(x)) obtained for the graphs with n = 6 (open circles) 
is the closest one to the Gaussian distribution (solid line) which is the prediction of RMT. 
Thus, the behavior of I n and P(ty(x)) for graphs with n > 6 suggest that the counterintuitive 
results obtained for the autocorrelation function c(x), where the best agreement with RMT 
was obtained for the graphs with six vertices, are connected with the localization effects. 

The above results suggest that in order to describe numerically the behavior of the auto- 
correlation function c(x) one can use the concept of coupled M-GOE matrices [42J. A coupled 
M-GOE matrix with M x M blocks can be constructed using a random GOE matrix. The 
values of diagonal blocks elements of the random block matrix are equal to the values of 
random GOE matrix elements, while the elements of off-diagonal blocks are equal to random 
GOE matrix elements multiplied by a coupling parameter < A < 1. By varying A between 
1 and one can change the coupling between the diagonal blocks of the matrix. In this way 
a transition between GOE case (A = 1) and M-GOE case (A = 0) can be achieved. In our 
numerical study of the parametric level dynamics of coupled M-GOE matrices we applied 
the Hamiltonian model ( fTTT) . but instead of Hi and if 2 to be GOE matrices we used two 
400 x 400 coupled M-GOE matrices. 

Let us consider the autocorrelation function c(x) for graphs with n > 6 vertices. Figure [3] 
shows that the parametric dynamics of graphs with n = 6 is closely described by RMT 
for GOE, while large discrepancies observed for graphs with n > 6 suggests that their 
properties can be rather described by coupled M-GOE matrices. In order to check this 
assumption we fitted the values of A and M parameters to describe properly the behavior 
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FIG. 5: The level velocity autocorrelation function c(x) for quantum graphs with: n = 6 (open 
circles), n = 7 (full squares), n = 10 (open triangles) and n = 30 vertices (full circles) compared 
to the results of RMT (solid line) and coupled M-GOE matrices: M = 3, A = 0.08 (dotted line); 
M = 12, A = 0.05 (dash-doted line) and M = 50, A = 0.03 (dashed line). 

of the autocorrelation function c(x) for graphs with n > 7. In Figure [5] the results obtained 
for quantum graphs with n = 7, 10 and 30 are compared with corresponding curves for 
coupled M-GOE matrices with M = 3, A = 0.08; M = 12, A = 0.05 and M = 50, A = 0.03, 
respectively. In all cases good qualitative agreement between the results for quantum graphs 
(symbols in Figure [5]) and the corresponding ones for coupled M-GOE matrices (lines in 
Figure [5]) is observed. 

We also calculated the integrated nearest neighbor spacing distribution (INNS) for quan- 
tum graphs and coupled M-GOE matrices (see Fig. [6]). We found that INNS calculated for 
graphs with the size between n = 4 and n = 30 (in Fig. [H] we show only results for n=6 
and n=30) and for coupled M-GOE matrices (in Fig. [H] we show results for matrices with 
M = 50 and A = 0.03) is in good agreement with the prediction of RMT for GOE. This 
result is very interesting, while INNS and the velocity distribution calculated for quantum 
graphs are in agreement with RMT for GOE, the autocorrelation function c(x) departures 



10 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

s 

FIG. 6: The integrated nearest neighbor spacing distribution (INNS) for quantum graphs with 
n = 6 (open circles), n = 30 (full circles) and coupled M-GOE matrices with M = 50 and A = 0.03 
(crosses) compared to the results of RMT for GOE (solid line). 



from RMT predictions for GOE and depends on the graph size n. It seems that c(x) is much 
more sensitive to non-universal features of the spectrum and wave functions of the graphs 



such as eg., scars 



39 



40J than the other considered statistics. 
Another important quantity connected with parametric correlations is the velocity au- 
tocorrelation function c(u,x). From the practical point of view the definition given by the 
formula (j7]) is not suitable because it contains delta-function 5(x). The presence of delta- 
function means that one should calculate the correlations of level velocities of energy levels 
separated in energies by exactly the value of u> from one another. In practice the number of 
such levels is close to zero. That is why we substituted infinitely narrow delta-function by 
the weight function f(x, 5) with the finite width S, which allows to define an autocorrelation 
function c$(u,x) that is more suitable from the practical point of view 

£(/(&(* + x) - fc(x) - u, y/2S) &&4%l& 



cx(u;,x) = — = . (14) 

kl 

In our calculations we used a Gaussian weight function f(x, 5) 

/M) = ^exp(-^). (15) 
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FIG. 7: The velocity autocorrelation function c(w,x) for quantum graphs with (a) n = 6, (b) 
n = 10 and (c) n = 30 compared to the results of RMT for GOE (solid lines) and coupled M-GOE 
matrices (broken lines). The autocorrelation function c(oj,x) is calculated for four values of the 
parameter to: uj = (circles), to = 0.1 (triangles), u> = 0.25 (diamonds), and to = 1 (squares). 

The autocorrelation function cs(uj, x) defined by the formula f|T4|) allows to calculate corre- 
lations of level velocities of energy levels separated mutually by uj ± 5. 

Figures [7^, [7b and show the autocorrelation function of level velocities cs(u, x) for 
quantum graphs with n = 6, n = 10 and n = 30 vertices, respectively. Calculations were 
performed for four values of the parameter u = 0, 0.1, 0.25 and 1.0. Similarly to the paper [5] 
the parameter 5 was chosen to be equal to 0.03. The numerical curves obtained for quantum 
graphs are compared with the theoretical ones for GOE in the case of graphs with n = 6 
(Fig. [TJi) and for coupled M-GOE matrices for graphs with n = 10 and n = 30 (Fig. [7b and 
Fig. [7b, respectively). In each case the statistical averaging was performed over 99 graphs 



12 



configurations. 

We found the autocorrelation function c$(uj,x) calculated for graphs with n = 6 vertices 
(Fig.[7Ji) to be in good agreement with the theoretical prediction for GOE (solid lines) for all 
the values of the uj parameter. For quantum graphs with n = 10 and n = 30 vertices (Fig. [7b 
and Fig. [7b) the numerical results are in agreement with the ones obtained for coupled M- 
GOE matrices (broken lines in Fig. [7b and Fig. [7b) for most values of the parameter x. This 
fact confirms that coupled M-GOE matrices can be successfully used for the description of 
the parametric dynamics of quantum graphs. However, for some values of x small deviations 
from coupled M-GOE matrices results are observed. In Fig. [TJi and Fig. [7b the biggest 
deviation is observed for uj = 0.25 in the vicinity of the maximum of the autocorrelation 
function cs(0. 25, x). 

Let us consider now the level curvatures «, the second derivative of the unfolded energies: 

S. - f|. (16) 
We rescaled the curvatures for quantum graphs using the variance <r v : 

m = ~~- (17) 

The distributions of level curvatures are presented in Fig. Eh, Fig. [8b and Fig. Eb for graphs 
with n = 6, n = 10 and n = 30, respectively. It is visible that for all of the cases there 
is no agreement of the obtained results with the GOE prediction (solid line). Only for the 
graphs with n = 6 vertices the tail of the distribution can be described correctly by the 
GOE prediction. For the other graphs the level curvature distributions departure from the 
GOE one. Therefore, in the case of graphs with n — 10 and n = 30 we compare the level 
curvature distribution with the results obtained for coupled M-GOE matrices (broken line). 
For the graphs with n — 10 vertices (Fig. [8b) small deviations from the results for coupled 
M-GOE matrices are observed only for small values of level curvatures. In the case of graphs 
with n = 30 (Fig. [Hb) the curvature distribution calculated for coupled M-GOE matrices 
(broken line) follows the results obtained for the graphs for all of the range of the parameter 

K. 

One of the characteristic features of level dynamics of quantum chaotic systems are 
avoided crossings which are often observed between the neighboring levels when an ex- 
ternal parameter X is changed (see Fig. [1] and Fig. [2]). We define the avoided crossing 
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FIG. 8: The parametric curvature distribution P(k) for quantum graphs with (a) n = 6 (open 
circles), (b) n = 10 (open triangles) and (c) graphs with n = 30 vertices (full circles). The 
numerical results are compared to the results of RMT for GOE (solid line) and for coupled M- 
GOE matrices (broken line). Results are presented in double logarithmic scale. 

distance C as a local minimum of the distance between two neighboring eigenvalues. The 
distribution of avoided crossings P(c) can be used to distinguish between chaotic systems 
which belong to different universality classes. The analytical formulas for the distributions 
of avoided crossings P(c) for GOE and GUE were derived by Zakrzewski and Kus [ljj]. For 
GOE-systems the distribution of avoided crossings P{c) is the following: 



To calculate the local minimum C of the avoided crossing distance we used the hyperbolic 




(18) 
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FIG. 9: The distribution of avoided crossings P(c) for quantum graphs with n = 6 (open circles), 
n = 30 (full circles) and coupled M-GOE matrices with M = 50 and A = 0.03 (crosses) compared 
to the results of RMT for GOE (solid line). 



formula 



Q: 



c 



\ 



(4 - d\f 



(19) 



•((/•J -\~ (l\ — ^<<2 j 

where tij, z = 1,2,3 are three consecutively calculated distances between the adjacent levels. 
Then avoided crossing distances C obtained in this way are rescaled using the formula: 

C = W (20) 

where (C) is the mean value of avoided crossing distance. 

Figure M shows the distribution of avoided crossings P(c) for quantum graphs with n = 6 
(open circles), n = 30 (full circles) and coupled M-GOE matrices with M = 50 blocks 
(crosses) compared to the results of RMT for GOE (solid line). The numerical results 
for graphs were obtained as a result of an averaging over 99 graph configurations. The 
distributions of avoided crossings P(c) obtained for graphs show an excess of small values 
compared to the theoretical prediction for GOE (solid line). However, it is worth noting 
that in the range of 0.8 < c < 1.5 the results obtained for the graphs with n = 6 vertices 
are slightly smaller than the RMT prediction for GOE. The results obtained for coupled 
M-GOE matrices show also an excess of small values comparable to the GOE prediction. 
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However, for c > 0.5 M-GOE results are close to the ones obtained for the RMT prediction 
for GOE as well as for the graphs with n = 30 vertices. 

In summary, we investigated numerically eigenvalue dynamics of fully connected irregu- 
lar quantum graphs with Neumann boundary conditions of different size n = 4 — 30. We 
showed that there are some spectral statistics such as the integrated nearest neighbor spac- 
ing distribution, the parametric velocity distribution P{y ) and the distribution of avoided 
crossings P(c) that show no or weak sensitivity to the system size. These statistics show 
GOE-like behavior. There are also some other statistics, e.g., the second order level velocity 
autocorrelation functions c(x) and c(u,x) and the parametric curvature distribution P(k) 
that for larger graphs (n > 6) show deviations from the predictions of RMT for GOE. In 
all these cases the obtained results are much better described by the model of coupled M- 
GOE matrices. The agreement of the obtained results with the coupled M-GOE matrices 
predictions supported by the results of the inverse participation ratio and the distributions 
P(^(x)) suggest the existence of non-ergodic structures of graphs' wave functions connected 
with localization which cause the quantum dynamical disintegration of large graphs into the 
union of weakly interacting smaller graphs. 

It is worth pointing out that the agreement of parametric spectral statistics of quantum 
graphs with the predictions of RMT could be restored by making quantum graphs "more 
ergodic". It could be done by randomizing the boundary conditions at the vertices. For 
example, one can replace vertex-scattering matrices by random matrices taken from Circular 
Orthogonal Ensemble (COE) 43]. Our preliminary results, which are not presented here, 
suggest that in this case parametric spectral statistics show good agreement with RMT 
predictions for GOE both for smaller and larger graphs. Moreover, this agreement becomes 
better for larger graphs. 
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